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Abstract. We present early time high-resolution (VLT/UVES) and late time low-resolution (VLT/FORS) optical spectra of the 
normal type la supernova, SN 2001el. The high-resolution spectra were obtained 9 and 2 days before (B-band) maximum light. 
This was in order to allow the detection of narrow hydrogen and/or helium emission lines from the circumstellar medium of 
, the supernova. No such lines were detected in our data. We therefore use these spectra together with photoionisation models to 

derive upper limits of 9 x 10~ 6 M Q yr~' and 5 x 10 5 M yr~' for the mass loss rate from the progenitor system of SN 2001el 
assuming velocities of 10 km s and 50 km s -1 , respectively, for a wind extending to outside at least a few xlO 15 cm away 
fT") , from the supernova explosion site. So far, these are the best Ha based upper limits obtained for a type la supernova, and exclude 

a symbiotic star in the upper mass loss rate regime (so called Mira type stars) from being the progenitor of SN 2001el. The 
low-resolution spectrum was obtained in the nebular phase of the supernova, ~400 days after the maximum light, to search 
for any hydrogen rich gas originating from the supernova progenitor system. However, we see no signs of Balmer lines in our 
spectrum. Therefore, we model the late time spectra to derive an upper limit of ~0.03 M Q for solar abundance material present 
at velocities lower than 1000 km s within the supernova explosion site. According to numerical simulations of Marietta et 
al. (2000) this is less than the expected mass lost by a subgiant, red giant or a main-sequence secondary star at a small binary 
f separation as a result of the SN explosion. Our data therefore exclude these scenarios as the progenitor of SN 2001el. Finally, 

JL , we discuss the origin of high velocity Ca II lines previously observed in a few type la supernovae before the maximum light. 

We see both the Ca II IR triplet and the H&K lines in our earliest (-9 days) spectrum at a very high velocity of up to ~34 000 
km s -1 . The spectrum also shows a flat-bottomed Si II '6150 A' feature similar to the one previously observed in SN 1990N 
(Leibundgut et al. 1991) at 14 days before maximum light. We compare these spectral features in SN 2001el to those observed 
in SN 1984A and SN 1990N at even higher velocities. 
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1. Introduction 

While the progenitor stars for a few core-collapse supernovae 
(SNe), viz. SN 1987A (White & Malin 1987; Walborn et al. 
1989), SN 1993J (Aldering et al. 1994; Maund et al. 2004) 
and SN 2003gd (Smartt et al. 2004), have been directly de- 
tected in pre-explosion images, and a general good picture 
has emerged from stellar evolutionary modelling (Smartt et al. 
2003; Eldridge & Tout 2004), similar constraints are still miss- 
ing for the progenitors of thermonuclear (type la) supernovae 
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(henceforth SNe la). SNe la are believed to arise from the ther- 
monuclear explosion of a carbon oxygen white dwarf (CO WD) 
in a binary system when its degenerate mass becomes close to 
the Chandrasekhar limit. Models predict that the companion 
star can either be another CO WD or a non-degenerate star. In 
the double degenerate scenario two low-mass CO WDs merge, 
whereas in the single degenerate scenario matter is accreted 
from the companion star until the CO WD reaches its limiting 
mass for carbon ignition (see Hillebrandt & Niemeyer 2000 
and Nomoto et al. 2000 for discussion and references). In nei- 
ther scenario would the progenitor system be luminous enough 
to be directly detectable at extragalactic distances. However, 
the recently reported (Ruiz-Lapuente et al. 2004) likely identi- 
fication of the subgiant/main-sequence companion star surviv- 
ing the explosion of Tycho Brahe's supernova (SN 1572) gives 
support for the single degenerate scenarios. 
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In all the single degenerate scenarios hydrogen and/or he- 
lium rich gas originating from the companion star could be 
present in the circumstellar medium (CSM) of the SN. This 
gas could be in the form of a stellar wind, an accretion disk, 
a filled Roche-lobe, a common envelope around the progenitor 
system, or a wind from the accreting WD (Hachisu et al. 1999a, 
1999b). This material could be detectable by making use of the 
SN as a probe of its own CSM. If the wind is dense enough and 
ionised by the radiation from the SN ejecta/wind interaction, 
narrow emission lines similar to those in narrow-line (type Iln) 
core-collapse SNe (and SN 1987 A at late times) could be visi- 
ble on top of the SN spectrum (Cumming et al. 1996; Lundqvist 
et al. 2006). In fact, SN 2002ic was recently observed to show 
spectral features similar to both SNe la and SNe Iln (Hamuy 
et al. 2003; Baron 2003; Deng et al. 2004; Kotak et al. 2004); 
it could belong to the so-called type 1.5 SNe (e.g., Chugai et 
al. 2004; Nomoto et al. 2004) which were predicted by Iben 
& Renzini (1983). The SN la explosion itself is also expected 
to cause a loss of material from a close, non-merging compan- 
ion (Fryxell & Arnett 1981; Taam & Fryxell 1984; Livne et al. 
1992; Marietta et al. 2000). This material is likely to be hydro- 
gen rich, and could reveal itself when the ejecta have become 
transparent. 

In addition, Ca II IR triplet (8498 A, 8542 A, 8662 A) and 
H&K (3968 A, 3934 A) lines have recently been observed at 
very high velocities compared to the SN photospheric lines 
in a few SNe la: SN 1994D (Hatano et al. 1999), SN 1999ee 
(Mazzali et al. 2005a), SN 2000cx (Thomas et al. 2004; Branch 
et al. 2004), SN 200 lei (Wang et al. 2003; Kasen et al. 2003), 
SN 2003du (Gerardy et al. 2004) and SN 2004dt (Wang et al. 
2004). While the origin of these high velocity (HV) lines is still 
unknown, one of the interpretations has been interaction of the 
SN ejecta with CSM material originating from the companion 
star (e.g., Gerardy et al. 2004). 

To shed light on the origin of SNe la we have been conduct- 
ing a Target of Opportunity (ToO) programme on the European 
Southern Observatory (ESO) Very Large Telescope (VLT). We 
have observed two SNe la so far, SNe 2000cx and 200 lei with 
the VLT (for preliminary analyses see Lundqvist et al. 2003; 
Lundqvist et al. 2005). The observations of SN 2000cx, and 
their modelling, are described in Lundqvist et al. (2006). Here 
we report our observations of SN 200 lei. 

Supernova 2001el was discovered on September 17.1 UT 
(Monard 2001) about 14" west and 20" north of the nucleus 
of an edge-on spiral galaxy NGC 1448. This nearby (V rec = 
1 164 km s _1 ; de Vaucouleur et al. 1991) spiral has hosted also 
two type II SNe (SN 1983S and SN 2003hn) during the last 
20 years. The magnitude of SN 2001el at the time of the dis- 
covery was about 14.5. It was classified as a SN la by us on 
September 21.3 UT (Sollerman et al. 2001). SN 2001el was 
observed to show a polarised spectrum before the maximum 
light (Wang et al. 2003) indicating some asymmetric structure 
during these initial phases of expansion. Later on, SN 200 lei 
has been well monitored both photometrically and spectro- 
scopically, and has been shown to be a normal SN la with 
Ami 5 (B)=1.13 (Krisciunas et al. 2003). 

In Sect. 2 the VLT observations and data reductions are 
described. In Sect. 3, we describe how we searched for CSM 



lines from the early time high-resolution and late time low- 
resolution spectra of SN 200 lei. No CSM lines were detected 
and thus upper limits were derived for the line fluxes. In Sect. 4, 
we use these flux limits together with models to estimate up- 
per limits for the amount of hydrogen rich material lost from 
the companion star either via a stellar wind or as a result of 
the SN explosion. In Sect. 5, we compare the profiles of the 
high velocity Ca II IR triplet and H&K lines as well as the 
flat-bottomed Si II '6150 A' feature in our early time spectrum 
to those observed in SNe 1984A and 1990N at a similar early 
epoch. Finally, in Sect. 6 we discuss our results and present the 
summary. 

2. VLT observations and data reductions 

2.1. High-resolution UVES spectroscopy before the 
maximum light 

Two days after the discovery IAU circular (Monard 2001) on 
September 21.3 (UT), SN 200 lei, still without a spectroscopic 
classification, was observed with the Ultraviolet and Visual 
Echelle Spectrograph (UVES) on the VLT. All the observations 
were obtained in service mode. We obtained two 2400 sec- 
ond exposures in each of two different dichroic settings, DIC1 
390+564 and DIC2 437+860 (Table 1). These settings together 
give a complete wavelength coverage between 3260 A and 
10600 A. The night was clear and the seeing during the obser- 
vations was around 0'.'8. The spectrophotometric standard star 
LTT 1020 (Hamuy et al. 1994) was also observed in both se- 
tups, after the supernova observations. All spectra were quickly 
reduced using the UVES pipeline. This allowed a classification 
of the supernova as a type la, well before the maximum light 
(Sollerman et al. 2001). 

The SN was observed again, 7 days later, on September 
28.3 (UT) to enable monitoring of any weak CSM lines po- 
tentially detectable from this promising early SN la event. 
Three exposures of 3000 seconds each were obtained using 
the dichroic setting DIC1 390+564, which covers the wave- 
length ranges 3260-4450 A and 4580-6680 A. The night 
was clear and the seeing was around l'.'l during the observa- 
tions. The spectrophotometric standard star HR 1996 (Tinney 
& Reid 1998) was observed immediately after the supernova 
observations. The 0V8 slit was used for all the observations on 
both epochs yielding a spectral resolution of R ~ 50 000 (~6 
km s _1 ) for the region around Ha. Our UVES observations of 
SN 2001el were thus obtained 8.7 and 1.7 days before the su- 
pernova B-band maximum light, which occurred on September 
30.0 (UT) (Krisciunas et al. 2003). 

2.2. UVES data reductions 

The spectral orders covering the lines of interest for this study 
were carefully reduced using a combination of tasks from the 
UVES pipeline 1 (Ballester et al. 2003) (version 2.0) and IRAF. 
All the 2D data were interactively reduced using the UVES- 
pipeline task REDUCE/SPAT and calibration frames obtained 

1 http://www.eso.org/observing/dfo/quality/ 
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Table 1. Log of VLT observations of SN 2001el. 



Instrument 


Date (UT) 


JD 


Epoch 


Exp. 


a.m." 


Seeing* 


set-up 


Slit width 






(2452000+) 


(days) 


(sec) 




(arcsec) 




(arcsec) 


UVES 


2001 Sep 21.22 


173.72 


-8.7 


2400 


1.30 


0.8 


390+564 c 


0.8 


UVES 


2001 SeD 21 25 


173.75 


-8.7 


2400 


1.19 


0.8 


390+564 


0.8 


UVES 


2001 Sep 21.28 


173.78 


-8.7 


2400 


1.12 


0.7 


437+860 d 


0.8 


UVES 


2001 Sep 21.31 


173.81 


-8.7 


2400 


1.08 


0.9 


437+860 


0.8 


UVES 


2001 Sep 28.22 


180.72 


-1.7 


3000 


1.20 


1.1 


390+564 


0.8 


UVES 


2001 Sep 28.26 


180.76 


-1.7 


3000 


1.11 


1.1 


390+564 


0.8 


UVES 


2001 Sep 28.29 


180.79 


-1.7 


3000 


1.07 


1.1 


390+564 


0.8 


FORS1 


2002 Nov 02.15 


580.65 


+398 


3000 


1.15 


1.2 


GRIS.300V 


1.3 



" The effective airmass of the observation 

b Seeing from the DIMM-monitor 

c Setting 390+564 covers wavelength ranges 3260-4450 A, 4580-6680 A. 

d Setting 437+860 covers wavelength range 3730-4990 A, 6600-10600 A. 



during daytime operation. This included bias subtraction, flat 
fielding the data in pixel space, automatically determining the 
Echelle order positions and slopes, and subtraction of the inter- 
order (median) background. The pipeline performed wave- 
length calibration for the 2D data by automatically identifying 
and fitting a large number (~ 1300 per set-up) of ThAr arc lines. 
For the wavelength calibration a single dispersion solution, ob- 
tained for the center coordinate of the object, is used for each 
pixel in the spatial direction. The pipeline then produces sepa- 
rate 2D frames for each Echelle order with constant wavelength 
binning on the x-axis, and the spatial extent of the data on the 
y-axis. 

An accurate correction for the background emission from 
the supernova host galaxy is vital when searching for traces of 
weak emission or absorption lines from the supernova CSM 
(see also Delia Valle et al. 1996). In fact, the widths and lumi- 
nosities of emission lines expected from the CSM of a SN la 
are similar to the ones observed from ordinary H II regions. The 
luminosity time evolution of such circumstellar emission lines 
can also be slow enough (in a timespan of a few days) to re- 
sult in confusion with the constant emission from a nearby H II 
region. We performed the background subtraction outside the 
UVES pipeline using the IRAF BACKGROUND task interactively. 
We found that the dispersion solution obtained for the center 
coordinate of the object was valid over the slit length (8-12") 
enabling accurate subtraction of the sky and host galaxy back- 
ground emission from the wavelength calibrated frames. 

The UVES pipeline uses a Gaussian approximation for the 
cross-order light profile during the optimal extraction. This is 
not valid when the signal to noise ratio (S/N) is high enough 
(Piskunov & Valenti 2002; Ballester et al. 2003). In fact, we 
found that the standard optimal extraction routines included, 
e.g., in the UVES pipeline and IRAF did not give satisfactory 
results for our Echelle spectra with high S/N. Therefore, in this 
study we decided to use simple average extraction for the SN 
and standard star spectra, and performed the cosmic ray rejec- 
tion when combining the ID spectra. 

To correct for the telluric absorption lines and any resid- 
ual instrumental features left after the flat fielding, the super- 
nova spectra were divided by the spectrum of the standard 
star observed during the same night as the supernova with the 



same instrument settings and scaled to match the effective air- 
mass of the science observation. Both the standard stars ob- 
served, LTT 1020 and HR 1996, show broad Ha absorption 
lines with their wings extending close to the region of inter- 
est when searching for the circumstellar Ha line from the su- 
pernova spectrum. We therefore normalised the wavelength re- 
gions of interest in the standard star spectra with a low order 
polynomial fit before the division of the supernova data. The 
LTT 1020 spectrum was also smoothed before the division in 
order to conserve the S/N of the science data. The corrected 
science exposures were then normalised by a low order polyno- 
mial fit, and average combined together (see Figs. 3-4). For this 
we used weighting by the inverse squares of the corresponding 
pixel to pixel noise levels. Any deviant pixels or cosmic rays 
(above 5cr level) were rejected and replaced by the average of 
the two unaffected spectra or in the case of only two available 
spectra by a fit to the continuum. The individual orders cover- 
ing the wavelengths around the Ca II (H&K and the IR triplet) 
and the Si II '6150 A' lines were carefully merged together in 
IRAF (see Fig. 8). 

The absolute flux calibration was performed for each wave- 
length region of interest separately. For this we used the well- 
sampled V-band light curve of SN 200 lei from Krisciunas et al. 
(2003) and a low-resolution spectrum from Wang et al. (2003). 
This spectrum was observed on September 26.3 (UT), i.e., 4.0 
days before the maximum light. To estimate the uncertainties 
due to the spectral evolution between the epochs of our UVES 
spectra and the spectrum of Wang et al. (2003) we used low- 
resolution spectra of SN 2002bo (Benetti et al. 2004) from 9.0, 
4.0 and 2.0 days before the maximum light. These are both nor- 
mal SNe la with exactly the same Ami 5 values (Krisciunas et 
al. 2003; Benetti et al. 2004). The low-resolution spectra were 
scaled in flux to match the observed V-band photometry at the 
dates of our UVES observations. From this comparison we esti- 
mate that our absolute flux calibration is accurate to better than 
+25%. 

The reduction procedure outlined above was designed to 
mainly study limited wavelength regions around narrow, and 
presumably weak line features. To obtain an overall spectrum 
of the SN we do not need such a high accuracy. For the first 
epoch (-8.7 day) overall spectrum we therefore simply used 
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Fig. 1. UVES spectrum of SN 2001el obtained on September 21.3 (UT), i.e., 8.7 days before the B-band maximum light. The 
spectrum has been median filtered and rebinned to have a resolution of about 50 km s _1 pixel -1 . The lines of interest for this 
study are indicated. The spectrum of SN 1990N observed (Leibundgut et al. 1991) at 14 days before the B-band maximum light 
is also shown for comparison. This spectrum has been redshifted to match the recession velocity of SN 2001el, scaled by x3.5, 
and shifted upward by 1.5xl0 -14 erg s" 1 cm^A^ 1 for clarity. The Ca II IR triplet and the H&K lines are present at a very high 
velocity (HV), and the characteristic Si II '6150 A' absorption shows a similar strongly blueshifted flat-bottomed profile in both 
SNe. However, we see no signs of narrow circumstellar hydrogen or helium lines nor HV H/3 in the SN 200 lei spectrum. The 
Ca II and Na I ISM lines are also marked in the SN 2001el spectrum (for details see Sollerman et al. 2005). A large number of 
telluric lines are present, especially at wavelengths longer than 9000 A. 



the standard UVES pipeline task REDUCE/UVES. Here the in- 
dividual orders were automatically extracted and merged to- 
gether by the pipeline. To enable successful automatic merg- 
ing, simple average extraction was used, and the individual or- 
ders were flat fielded in ID. The relative flux calibration was 
carried out in IRAF separately for the spectra corresponding 
to each spectrograph set-up (see Table 1). For this we used 
the standard star (LTT 1020) spectra observed together with 
the SN. To remove the large number of cosmic rays we me- 
dian filtered each individual spectrum using a filter window of 
50 pixels (~50 km s _1 ). The overlapping spectra were average 
combined together, and rebinned to yield a pixel scale of ~50 
km s _1 pixel -1 . The resulting spectrum is shown in Fig. 1 where 



the levels of the wavelength regions covered by the different 
set-ups have been adjusted by eye to match the overall shape of 
the -4 day low-resolution spectrum of Wang et al. (2003). The 
overall level of the spectrum has been scaled in flux to match 
the observed V-band photometry of the SN. 

2.3. Low-resolution FORS spectroscopy in the nebular 
phase 

We have also obtained a very late low-resolution optical spec- 
trum of SN 200 1 el. On November 2.2 (UT) 2002, 398 days past 
maximum light, we used FORS 1 on the VLT to obtain a 3000 
second exposure (see Table 1 for details). We used the 300V 
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Fig. 2. Ca II H (,13968, dashed) and K (,13934, solid) inter- 
stellar absorption lines towards SN 200 lei within its host 
galaxy NGC 1448. The velocity scale is heliocentric. See 
Sollerman et al. (2005) for details about the absorption com- 
ponents. 

grism together with order-sorting filter GG375 and a 1'.'3 wide 
slit yielding a spectral resolution of R ~ 440 (~700 km s _1 ) for 
the observation. 

The spectrum (see Fig. 6) was reduced in a standard way in 
IRAF, including bias subtraction, flat fielding, and wavelength 
calibration using spectra of a Helium-Argon lamp. Flux cal- 
ibration was done relative to the spectrophotometric standard 
star LTT 2415 (Hamuy et al. 1994). We have also compared 
the synthetic V band magnitude obtained from this spectrum 
to late time light curves of SNe la (e.g., Sollerman et al. 2004, 
their Fig. 11) and conclude that the absolute fluxing of the spec- 
trum is accurate to +20%. 

3. Searching for traces of circumstellar lines 

3.1. The expected wavelength range for the CSM lines 

Superposed on the supernova spectra (Fig. 1) we also detect in- 
terstellar Ca II H&K absorption both from the Galaxy and from 
NGC 1448 which hosted the supernova. The Ca II H&K ab- 
sorption lines arising within the host galaxy are plotted in Fig. 2 
as observed on Sept. 28. Several absorption components are 
present in the spectrum (see detailed analysis in Sollerman et 
al. 2005). The overall range of velocities for the absorbing gas 
along the line-of-sight to the supernova is V a t, s ~ 1130-1230 
km s _1 . The recession velocity of the H II regions close to the 
supernova included in the UVES slit can be used as another in- 
dicator of the supernova recession velocity. We measured Vh ii 
= (11 86 ± 45) km s _1 for the Ha emitting gas next to the SN po- 
sition in the UVES spectrum at Sept. 28. As NGC 1448 is ori- 
ented almost edge on (z = 88°) the galaxy rotation may result in 
a substantial velocity for the supernova w.r.t. to the host galaxy 
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Table 2. Comparison of limiting fluxes from two different 
methods. The center of the Gaussian noise distribution, fj., 
and its width, cr, yield 3cr upper limit for the emission line 
profile peak flux density, f(3cr)/. These can be compared to 
the f(3cr) /7 values derived simply from the measured pixel-to- 
pixel noise in adequately rebinned spectra. The parameters of 
the flux density distributions (fi and <x) and the upper limits 
are given in units of 10~ 17 erg s _1 cirT 2 A -1 . 



Line 


Epoch 


FWHM 




<r 


1(30-), 


f(3<r) 






(kms- 1 } 


(io- 


17 

erg 


-1 -2 

s cm 


A" 1 ) 


Ha ,16563" 


1 


37 


0.17 


2.9 


8.9 


11 


Ha ,16563 


1 


62 


-0.67 


2.6 


7.1 


8.1 


Ha ,16563 


2 


37 


1.1 


3.6 


12 


11 


Ha ,16563 


2 


62 


1.2 


2.4 


8.4 


8.8 


He 1 ,15876* 


1 


21 


0.34 


7.7 


23 


29 


He 1 ,15876 


1 


53 


1.7 


5.4 


18 


15 


He II /I4686 c 


1 


21 


0.06 


3.9 


12 


14 


He II ,14686 


1 


53 


0.07 


2.2 


6.7 


7.2 


He II ,14686 


2 


21 


1.3 


7.7 


24 


23 


He II ,14686 


2 


53 


1.2 


5.0 


16 


16 



" The Ha upper limits were derived for the wavelength range 
6580 A- 6597 A. 

b The He I upper limits were derived for the wavelength range 
5891 A- 5910 A. 

c The He II upper limits were derived for the wavelength range 
4695 A- 47 13 A. 



recession velocity. A rotation velocity of V 10t ~193 km s 
has been measured for NGC 1448 (Mathewson & Ford 1996). 
We therefore carry-out the search for the narrow circumstellar 
emission lines at the range of velocities V lec = (1180 + 190) 
km s -1 . This is close to the recession velocity of NGC 1448 of 
1164 + 5kms~' compiled by de Vaucouleur et al. (1991). 

3.2. Early time spectra 

The overall UVES spectrum obtained on Sept. 21.3 (UT), i.e., 
8.7 days before the B-band maximum light is shown in Fig.^ 
The spectrum of SN 1990N observed at a very early epoch of 
~14 days before the B-band maximum light (Leibundgut et al. 
1991) is also shown for comparison. Although SN 1990N was 
at a much earlier phase than SN 2001el at the time of the ob- 
servation these spectra appear remarkably similar. Both SNe 
show strong Ca II IR triplet and H&K lines at similar high ve- 
locities. The characteristic Si II '6150 A' absorption appears 
flat-bottomed in both the spectra. To our knowledge such flat- 
bottomed Si II lines have not been seen in any other SNe la. 
The origin of these lines will be discussed in Sect. 5. 

The strongest lines expected to form in a hydrogen or 
helium-rich companion wind, and which are covered by our 
spectral range, i.e., Ha, or He I ,15876 and He II ,14686 (cf. 
Lundqvist et al. 2006), were searched for in the UVES spectra 
at the two epochs of observation (see Figs. 3-4). For the pre- 
dicted temperatures of -20 000 - 40 000 K (Lundqvist et al. 
2006) thermal broadening dominates the line widths for wind 
velocities up to ~40 km s _1 for hydrogen and ~20 km s _1 for 
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Fig. 3. Normalised and rebinned (~4 km s _1 pixel 1 ) UVES spectra in the expected spectral region around Ha for SN 2001el on 
two epochs, September 21.3 and 28.3 (UT) 2001, i.e., 8.7 and 1.7 days before the SN maximum light, respectively. The expected 
wavelength range of Ho- is marked with a horizontal dashed line, and the upper spectrum has been shifted vertically for clarity. 
No significant emission or absorption lines are visible. 
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Fig. 4. Normalised and rebinned (~4 km s 1 pixel ') UVES spectra in the expected spectral regions around He II (4686 A) and 
He I (5876 A) for SN 2001el on September 21.3 (UT) 2001, i.e., 8.7 days before the maximum light. The expected wavelength 
range of the CSM lines is marked with a horizontal dashed line. No significant emission or absorption lines are visible. 



helium. Here we will assume an average temperature of 28 000 
K (for details see Lundqvist et al. 2006) and a wind velocity, 
V„i n d, between 10 km s _1 and 50 km s _1 , and therefore adopt 
observed line widths of 37-62 km s _1 and 21-53 km s _1 for 
hydrogen and helium lines, respectively. 

No narrow lines with the expected line widths were de- 
tected in our spectra and therefore only upper limits for the 
emission line fluxes were obtained. To estimate detection lim- 
its for the narrow CSM emission lines we used the IRAF 
FITPROFS task to perform least squares fitting of Gaussian 
profiles (one at a time) to the normalised and unbinned spec- 
tra. The free parameters were the amplitude, the width (sigma) 



and the wavelength position of the Gaussian. As the spectra 
were already normalised, we did not fit the background. The 
lcr width of the Gaussian profile was fixed according to the 
relevant FWHM, and the line center was moved through the 
spectrum in small (FWHM/2) steps, i.e., only the amplitude 
of the profile was kept as a free parameter. To obtain statisti- 
cally meaningful upper limits lines were fitted to a +2 x 190 
km s _1 region around the most probable location of the CSM 
line with similar noise characteristics. This procedure yields 
a list of potential emission and absorption lines through the 
spectrum. In order to derive reliable flux limits from the data, 
the histogram distributions of the peak flux densities of the line 
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Fig. 5. Histograms of the line profile peak flux density distributions in the the first epoch UVES spectrum at the wavelength range 
(+2 x 190 km s ) around the expected position of Ha and assuming a line profile with FWHM of 37 km s (upper-left) and 
62 km s _1 (upper-right). The histograms have been fitted with Gaussians yielding reduced x 2 $ of 3.2 and 1.8, respectively. The 
centers (/i) and widths (cr) of the distributions are given in units of 10~ 17 erg s _1 cirT 2 A -1 . Lower-left: The unbinned spectrum 
(0.77 km s _1 pixel -1 ) showing the wavelength region around the expected position of Ha in our first epoch spectrum. A Gaussian 
line profile with FWHM of 37 km s _I and peak equal to the 3cr level is also shown. Lower-right: The same as on the left except 
that now the spectrum has been rebinned to ~18 km s _1 pixeL 1 resolution. 



profiles were inspected. These histograms have centers close 
to zero flux density and can be fitted with Gaussian profiles 
(see Fig. 5 (upper panel)) yielding reduced x 2 s of around 2- 
3 as expected from pure random noise. The centers (/i) and 
lcr widths of these distributions together with the derived 3<x 
emission line upper limits (f(3cr) = fi + 3 x cr) are listed in 
Table 2. For comparison we also show 3cr upper limits derived 
simply from the measured pixel-to-pixel standard deviations in 
the spectra rebinned to have a Nyquist sampling, i.e., one pixel 
equal to FWHM/2 for each expected line width (see above). 
This simple approach yields 3<x limits consistent with the more 
mathematically robust method described above. 

In Fig. [5] (lower panel) we show the wavelength region 
around the expected position of the Ha line in our first epoch 
spectrum together with a Gaussian profile with FWHM = 37 
km s and its peak at the 3cr level as obtained using the method 
described above. For lines of 20 - 60 km s the reduced UVES 
data with a 0.77 km s _1 pixel size are strongly over-sampled. 



Therefore, a comparison of the upper limits derived above with 
sufficiently rebinned data is more meaningful. In Fig.[5](lower 
panel) we show both the unbinned (0.77 km s pixeL 1 ) spec- 
trum and a spectrum rebinned to 18 km s _1 pixeL 1 . We note 
that there are no features with fluxes at the level or higher than 
the estimated 3cr level present in the rebinned data at the ex- 
pected wavelength range of Ha. The derived 3<x upper limits 
for the hydrogen and helium lines at the two different epochs 
are given in Table 3. For the first epoch observation we also 
give Ha upper limits for higher wind velocities of 100 km s _1 
and 200 km s . These limits are simply based on the pixel-to- 
pixel standard deviations in the rebinned (1 pixel = FWHM/2) 
spectra between 6570 A and 6608 A. 

3.3. Late time spectra 

In Fig.|5]we show the observed nebular spectrum of SN 2001el 
at 398 days after maximum light. The SN la spectra at these 
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Table 3. 3cr upper limits for circumstellar emission line 
fluxes (in units of 1(T 17 erg s -1 cm -2 ) assuming an average 
temperature of T = 28 000 K (see Lundqvist et al. 2006 for 
details). The upper limits are shown for the wind velocities 
ranging between 10 km s -1 and 50 km s _1 , i.e., the resulting 
line widths of FWHM = 37-62 km s~' and 21-53 km s~' for 
hydrogen and helium lines, respectively. For the first epoch 
observation the Ha upper limit is also given for the wind ve- 
locities of 100 km s" 1 and 200 km s" 1 . 



Line 


*wind 


FWHM 


Epoch 1 


Epoch 2 




(km s" 1 ) 


(km s" 1 ) 


(10- 17 erg 


; s" 1 cm -2 ) 


Ha ,16563 


10 


37 


<7.7 


<10 


Ha ,16563 


50 


62 


<10 


<12 


Ha ,16563 


100 


106 


<15 a 




Ha ,16563 


200 


203 


<20 fl 




He 1 ,15876 


10 


21 


<10 


_b 


He 1 ,15876 


50 


53 


<20 


_h 


He II ,14686 


10 


21 


<4.2 r 


<8.4 


He II ,14686 


50 


53 


<5.9 r 


<14 



The upper limits for V willd of 100 km s 1 and 200 km s 1 are based 
on the pixel to pixel standard deviations in the rebinned (1 pixel = 
FWHM/2) spectra between 6570 A and 6608 A. 
b Due to unsuccessful telluric correction the He I upper limits 
for the 2nd epoch are significantly higher than for the 1st epoch 
observation and therefore are not meaningful to present here. 
c The on-source exposure time for the wavelength region around 
He II /t4686 was twice as long as for the regions around the other 
lines observed in the first epoch since this region was covered by 
both the 437 and 564 settings of the UVES instrument. 



epochs are dominated by Fe emission. The features between 
4000 A and 5500 A are blends of Fe II and Fe III emission, 
while the features in the range 7000 A to 7600 A are mainly 
due to Fe II and Ni II. By comparing the observations to de- 
tailed late time modelling we are able to put a conservative 
upper limit on the amount of solar abundance material evapo- 
rated from a binary companion. Following the simulations of 
Marietta et al. (2000) we have searched in our spectra for signs 
of relatively narrow (FWHM < 1000 km s 1 ) Ha emission 
at the range of wavelengths corresponding to the recession ve- 
locity found in Sect. 3.1. No such emission is apparent in the 
spectrum. Our model calculations are described in Sect. 4.3. 

4. Constraints on the progenitor system of 
SN 2001 el 

4.1. Extinction and distance of SN 2001 el 

In order to estimate the upper limit for the mass loss rate, M, 
from the progenitor system of SN 200 lei we need to know the 
upper limits for the line luminosities, i.e., we need to know 
the extinction and distance towards the supernova. Krisciunas 
et al. (2003) have measured a colour excess of E(B - V) = 
0.253 + 0.063 for SN 2001el (see also Sollerman et al. 2005) 
using the Lira relation (Phillips et al. 1999). Here we adopt 
E(B -V) = 0.25 + 0.06, i.e., A v = 0.78+ 0.19 (assuming R v = 
3.1) for the extinction towards the supernova. 



Probably the most accurate method for deriving distances 
for normal SNe la is by making use of their photometry. 
Krisciunas et al. (2003) observed m(H) 14 = 13.18 + 0.05 for 
the SN at +14 days from the B-band maximum light. Meikle 
(2000) estimated M(H) 14 = -18.22+0.05 for "IR-normal" SNe 
la. Assuming the Galactic extinction law of Rieke & Lebofsky 
(1985) we have an H-band extinction of A H = 0.14 ± 0.03 to- 
wards the supernova. This results in a distance of 17.9+0.7 Mpc 
for the supernova. We note that Krisciunas et al. (2003) ob- 
tained exactly the same distance (17.9+0.8 Mpc) for SN 2001el 
using their BVI light curves and the method of Phillips et al. 
(1999). As the two distance estimates obtained by making use 
of the supernova photometry are identical we adopt 17.9+0.8 
Mpc for the distance of SN 200 lei. 

4.2. Upper limit for the progenitor mass loss rate 

We used our estimated upper limits for the narrow Ha line 
luminosity together with the photoionisation models from 
Lundqvist et al. (2006) to obtain upper limits for the mass 
loss rate from the supernova progenitor system. In these cal- 
culations wind velocities of 10 km s -1 and 50 km s" 1 and so- 
lar abundances were assumed. The model in Lundqvist et al. 
(2006) is an updated version of that described in Cumming et 
al. (1996). We still assume spherical symmetry and power-law 
density distributions of the SN ejecta and the circumstellar gas, 
where the ejecta profile is assumed to follow p e j e cta K r~ 7 , and 
the wind has the standard p w ; n d oc r 2 distribution. The shock 
structure, and blast wave expansion, can then be described by 
the similarity solutions of Chevalier (1982). The shocked gas 
produces hard radiation that ionises the unshocked wind and 
gives rise to narrow circumstellar lines. A difference compared 
to Cumming et al. (1996) is that we now assume much higher, 
and therefore more realistic, velocities of the ejecta (see also 
Sect. 6). The velocity of the fast circumstellar shock is in our 
present models V s ~ 4.45 X 10 4 km s -1 at day 1, which is 
about twice as fast as in Cumming et al. (1996). The use of 
such a high shock velocity is supported by the recent observa- 
tions of high velocity Ca II lines in early SN la spectra (see 
Section 5 and references therein) with their maximum veloci- 
ties up to ~35 000 km s _1 , especially since the circumstellar 
shock can be expected to advance ~ 1.2- 1.3 times (Chevalier 
1982) faster than the fastest unshocked ejecta that gives rise 
to the Ca II lines. In addition, the maximum velocity of un- 
shocked ejecta decreases with time, so the bluest absorption 
observed in the ejecta 1-2 weeks before the maximum light (see 
Table 3) can only give a lower limit to the maximum velocity 
of unshocked ejecta at even earlier times. A dense wind will 
therefore be overtaken much faster. Compared to Cumming et 
al. (1996) we have also calculated the free-free emission from 
the shocked ejecta more accurately (cf. Lundqvist et al. 2006). 

Adopting an extinction of Av = 0.78+ 0.19, a distance 
of D = 17.9 + 0.8 Mpc, and V wind = 10 km s 1 for the 
SN 200 lei progenitor system, we obtain 3cr upper limits for 
the Ha narrow line-luminosity of Lu a < 5.3 x 10 36 erg s 1 
and L Ha < 6.9 x 10 36 erg s" 1 for the first (-8.7 days) and the 
second (-1.7 days) epoch, respectively. If a higher wind veloc- 
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Fig. 6. FORS1 spectrum of SN 2001el obtained on November 2.2 (UT) 2002, i.e., 398 days past the B-band maximum light. 
The inset shows an enlargement of the region marked by dashed lines, focusing on wavelengths around Ha for the redshift of the 
SN. The SN emission at these wavelengths is dominated by a blend of Fe II emission lines (e.g., 6561 A). The dashed and dotted 
spectra in the inset are models calculated at 380 days for the SN, assuming 0.01 M and 0.05 M of solar abundance material, 
respectively, concentrated to +1000 km s _1 from the rest velocity of the supernova. This is to simulate Ha emission from gas 
removed from a possible companion star according to the models of Marietta et al. (2000). There is no sign of such Ha emission 
in SN 2001el at these wavelengths, which places a limit on hydrogen rich material of ~ 0.03 M . The model used to calculate 
the late emission is described in Sect. 4.3. 



ity of 50 km s is assumed, the luminosity upper limits are 
L Ha < 6.9 x 10 36 erg s _1 and L Ha < 8.3 x 10 36 erg s" 1 for 
the two epochs, respectively. In Fig. these upper limits are 
plotted together with modelled line luminosities for six differ- 
ent progenitor mass loss rates. Here we have assumed a rise 
time of 19.5 + 0.2 days for SN 2001el which is typical for 
normal (Ami 5 (B) =1.1) SNe la (Riess et al. 1999). Mass loss 
rate and luminosity are close to linearly proportional to each 
other in the logarithmic space (see Lundqvist et al. 2006 for 
details). Therefore the uncertainties in the rise time, distance 
and extinction have only a small effect in the derived mass loss 
rates. Assuming a wind velocity of 10 km s _1 our observed 
line luminosity limits therefore correspond to mass loss rates 
of about < 8.9 x 10~ 6 M Q yr 1 and < 1.2 x 10~ 5 M yr -1 , re- 
spectively. Somewhat higher upper limits of < 1.5 x 10~ 5 M 
yr -1 have previously been obtained for SNe 1994D (Cumming 
et al. 1996) and 2000cx (Lundqvist et al. 2006) also assuming 
Vwind = 10 km s _1 . However, if V w i„d = 50 km s _1 is as- 
sumed instead, the upper limits increase to about < 4.9 x 10~ 5 
M yr -1 and < 6.9 x 10~ 5 M yr -1 for the two epochs, re- 
spectively. The mass loss rates observed for symbiotic stellar 
systems are in the range between 10~ 8 M yr~' and 5 x 10~ 5 
M yr 1 (Seaquist et al. 1993) assuming V w i n d = 30 km s 1 . If 



adopting V w i n d = 10 km s or 50 km s _1 instead, the high- 
est mass loss rates observed for so called Mira type systems 
with a red giant secondary star would be ~2 x 10~ 5 M yr -1 
or ~8 X 10~ 5 M yr -1 , respectively. Therefore, our results are 
not consistent with a symbiotic stellar system in the upper mass 
loss rate regime being the progenitor of SN 200 lei. (See also 
Sect. 6.) 

For He I and He II the 3cr line luminosity upper limits are 
^He ~ 7 x 10 36 erg s _1 and < 3 x 10 36 erg s _1 , respectively, for 
the first epoch observation assuming a 10 km s _I wind. These 
are well above the line luminosities predicted by the photoion- 
isation models of Lundqvist et al. (2006). However, future ob- 
servations covering the He I /H0830 might give more meaning- 
ful results for a helium rich wind. 

4.3. Upper limit for the gas lost from the companion 
due to the SN explosion 

Marietta et al. (2000) presented two-dimensional simulations 
of the impact of a SN la explosion on a companion star. They 
studied three different scenarios in which they assumed the 
companion star to be a main-sequence star, a sub-giant or a 
red giant star. They found that the amount of of hydrogen rich 
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Fig. 7. The modelled evolution of the Ha narrow emission line luminosities are shown for the wind velocity of 10 km s (left) 
and 50 km s _I (right) as a function of the time relative to the B-band maximum light (assuming a rise time of 19.5 days for 
SN 2001el). The three different curves (dashed lines) correspond to mass loss rates of 5.6 x 10~ 6 M yr _1 , 1.0 x 10~ 5 M yr 
and 1.8 x 10~ 5 M Q yr" 1 , respectively in the LH figure, and 2.8 x 10~ 5 M yr"\ 5.0 x 10~ 5 M yr _1 and 9.0 x 10~ 5 M yr _1 in 
the RH figure. The 3<x upper limits for the Ha line luminosity at the two epochs of observation, -8.7 and -1 .7 days, are marked 
with a circle and triangle, respectively. 



material lost varied between ~0.002 M and ~0.5 M depend- 
ing on the companion star. They also found that most of the 
hydrogen-rich gas received velocities below ~1000 km s _1 in 
all three scenarios. In order to search for traces of hydrogen 
in the SN la ejecta as revealed by late time emission we have 
modelled nebular spectra using the code originally described 
in Kozma & Fransson (1998), and updated as described in 
Sollerman et al. (2004) and Kozma et al. (2005). The calcu- 
lations are based on the explosion model W7 (Nomoto et al 
1984; Thielemann et al. 1986) where we have artificially in- 
cluded varying amounts of solar abundance material. 

Based on the calculations by Marietta et al. (2000) we 
fill the inner region, out to 1000 km s _1 , in the W7 model 
with varying amounts of solar abundance material. Four time 
dependent calculations were performed including 0.01, 0.05, 
0.10, and 0.50 M of solar abundance material, (for further 
details see Lundqvist et al. 2006). As before an extinction of 
Ay = 0.78 and a distance of D = 17.9 Mpc were adopted for 
SN 2001el. Marietta et al. (2000) also studied the asymmetry 
of the ejecta due to the impact. The asymmetry effects were not 
taken into account in our one-dimensional spherically symmet- 
ric models. However, the asymmetries will mainly change the 
line profiles and will not affect our conclusions on the hydrogen 
emission. 

In Fig. [6] we show, together with the observations of 
SN 2001el, two modelled spectra at 380 days past explosion 
containing 0.01 M and 0.05 M of solar abundance material. 
We see that even in the model containing 0.05 M of hydrogen 
rich material, significant Ho- emission is present. This is clearly 
not seen in the observed spectrum. We therefore conclude that 
~0.03 M of solar abundance material would have been de- 



tected in our late time spectrum via relatively narrow (FWHM 
< 1000 km s~') Ha emission. Even such small amounts of 
solar abundance material show up in the nebular spectrum be- 
cause the optical depth to gamma-rays in the central, hydrogen 
rich region is sufficiently high to capture enough gamma-rays 
to power the radiation from hydrogen. This is in contrast to the 
Fe/Ni rich lower density regions that dominate the emission in 
the spectrum in Fig.[6]which are predominantly powered by lo- 
cal positron deposition. The mixing of newly synthesised ma- 
terial into the central regions would increase the deposition of 
positrons in the hydrogen rich material. Therefore, our derived 
limit on the mass of hydrogen gas from a companion in the 
central regions of the supernova ejecta should provide a con- 
servative upper limit. 



Marietta et al. (2000) found that the amount of stripped and 
evaporated gas from the companion depends on the nature of 
the star. The smallest amount of gas lost was for the main- 
sequence companions at the largest binary separations of six 
or more stellar radii. They lose < 0.02 M of solar abundance 
material due to the impact. However, the mass lost from the 
main-sequence companions at a smaller binary separation of 
~3 stellar radii, subgiant and red giant stars ranges from 0.15 
M to ~ 0.5 M , respectively. Most of this gas has a velocity of 
< 1000 km s _1 as the gas is evaporated from the companion. 
Scenarios with such high masses are clearly not favoured by 
our results. We emphasise that the amount of gas stripped off 
the companion at a high velocity is much lower than the evap- 
orated gas. High-velocity gas may only reveal itself in early 
spectra (see below). 
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5. High-velocity features in SNe 2001 el, 1990N, 
and 1984A 

Early spectroscopic observations of SNe la covering wave- 
lengths longer than ~ 7500 A are still scarce. Such observations 
have been reported to show the Ca II IR triplet lines (AA 8498, 
8542, 8662) at a high velocity (HV) compared to the SN photo- 
spheric velocities for a few SNe la. This includes SN 1994D at 
-8 days (Hatano et al. 1999), SN 1999ee at -7 days (Mazzali 
et al. 2005a), SN 2000cx at +2 days (Thomas et al. 2004), 
SN 2001el at -4 days (Wang et al. 2003; Kasen et al. 2003), 
SN 2003du at -5 days (Gerardy et al. 2004), and SN 2004dt at 
-7 days (Wang et al. 2004). 

Kasen et al. (2003) investigated a plausible geometry for 
the material giving rise to the observed HV Ca II IR triplet fea- 
ture in SN 200 lei and its spectropolarimetric properties (Wang 
et al. 2003). They favour a clumpy shell geometrically detached 
from the photospheric material, and suggest that the HV mate- 
rial could be Ca rich gas originating in either nuclear burning 
in the WD or from a swept up accretion disk surrounding the 
WD. 

Thomas et al. (2004) modelled both the HV IR triplet, as 
well as the HV H&K feature seen in SN 2000cx at +2 days. 
They found that 3D modelling was needed to produce satisfac- 
tory spectral fits for both the features. They investigated a sce- 
nario where the HV material would be hydrogen or helium rich 
gas in the SN CSM stripped from a non-degenerate companion 
star of the progenitor WD as suggested by Marietta et al. (2000; 
see also Sect. 4.3). The models of Thomas et al. (2004) require 
at least -10 -3 M of hydrogen or helium rich material in the 
HV shell, which they find a factor of ten higher than allowed 
by the simulations by Marietta et al. (2000) at such a high ve- 
locity. Due to the lack of a plausible donor mechanism to the 
SN CSM, Thomas et al. (2004) conclude that the HV material 
is more likely to have its origin in the SN explosion itself. 

Gerardy et al. (2004) further investigated the scenario of 
solar abundance material in the SN CSM. They carried out 
modelling where fast SN ejecta collide with stationary or slow 
moving CSM material close to the SN, setting up a forward 
and reverse shock structure. They assumed that the swept up 
CSM matter and shocked ejecta matter are mixed in the inter- 
action forming a dense shell of gas with a suitable Ca abun- 
dance (solar abundance) to give rise to the observed HV lines. 
In their modelling three different scenarios produced a shell of 
0.02 M at 20 days post-explosion which is capable of giving 
rise to the observed HV Ca II IR triplet feature. We note that 
the progenitor mass loss rate required by their first scenario is 
a factor of ~20 higher than our upper limit for the progenitor 
of SN 2001el (assuming a 10 km s _1 wind). This scenario as 
well as a scenario with a constant density CSM environment 
were also found unrealistic by Gerardy el al. (2004) as their 
energy generation rates would clearly dominate the bolomet- 
ric light curve of the SN. For SN 2003du Gerardy el al. (2004) 
favour a scenario where the matter involved in the interaction 
originates much closer to the progenitor system, in an accretion 
disk, Roche-lobe, or in a common envelope. 

More recently, Mazzali et al. (2005a) found that the ob- 
served HV Ca II IR triplet and Si II features in SN 1999ee could 



be explained by an increase of the mass at the highest ejecta ve- 
locities. According to the authors such a density increase could 
either be due to asymmetries in the explosion, and/or interac- 
tion of the ejecta with hydrogen rich material in the CSM of 
the SN. They found that about 4x10 3 M Q of hydrogen rich 
material in the CSM could explain the HV features observed in 
SN 1999ee. 

5.1. High-velocity Ca II 

In Fig. 8 we show Ca II profiles of SN 2001el together with the 
profiles of two other SNe la, SN 1984A (Wegner & McMahan 
1987) and SN 1990N (Leibundgut et al. 1991) both showing 
very high velocities at a similar early epoch. For these fig- 
ures the spectra were normalised by a first order polynomial 
fit around the profile of interest, and converted to the veloci- 
ties w.r.t. to the rest wavelengths assuming V rec = 1 180 km s _1 , 
1010 km s 1 and -261 km s 1 , for SNe 200 lei, 1990N, and 
1984A, respectively. To illustrate the maximum extent of the 
absorbing gas (in velocity space) we used the rest wavelengths 
of the bluest components of the Ca II triplet/doublet lines viz. 
8498 A and 3934 A. However, due to limited wavelength cov- 
erage in the SN 1990N spectrum the normalisation of its Ca II 
IR triplet profile is unreliable in the red part. Therefore, in the 
discussion below we will only make use of the blue part of this 
profile. In Table 4, we list the maximum velocities as indicated 
by the blue edges of the normalised HV profiles. 

Our first epoch (-9 days) observation of SN 200 lei consti- 
tutes to our knowledge the earliest reported coverage 2 of the 
HV Ca II IR triplet in a SN la (in addition to SN 1990N). 
Our -9 day observation covers also the Ca II H&K lines en- 
abling comparisons between the velocity profiles of these two 
features. In Fig. 8a we show profiles of the HV Ca II triplet and 
doublet lines as observed at -9 days. We illustrate the spectral 
evolution of the Ca II IR triplet between -9 and -4 days in 
Fig. 8b using the low-resolution spectrum obtained by Wang 
et al. (2003) for the latter epoch. At both epochs, the IR triplet 
profile of SN 200 lei shows a double dipped structure whereas 
the H&K doublet (3934, 3968 A) feature (-2600 km s _1 sepa- 
ration) does not show evidence for two components. The blue- 
most minimum in the IR triplet profile has been identified to be 
due to the blend of the /18542 and A8498 components (-1500 
km s _1 separation) whereas the redmost minimum is produced 
by the /I8662 component (Kasen et al. 2003). 

The velocities of the centers of the two Ca II profiles at -9 
days are quite similar. However, the maximum velocity, as in- 
dicated by the blue edge of the profile, is higher for the Ca II 
H&K lines (-34 000 km s 1 ) than for the IR triplet (-28 000 
km s~'). Similar behaviour was also shown by SN 1990N at 
-14 days with maximum velocities of 36 000 km s _1 and 
-28 000 km s _1 for the Ca II H&K and the IR triplet profiles, 
respectively. This behaviour can be due to blending of the Ca II 
H&K profile with other lines as already pointed out by Kasen 
et al. (2003) for SN 2001el, but may also just reflect a higher 

2 While this paper was being refereed other studies presenting 
data at even earlier epochs have appeared (see Mazzali et al. 2005b; 
Quimby et al. 2005). 
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absorption in the H&K lines. This can be expected since the 
H&K lines arise from a lower level which may well have a 
higher population than the level from which the Ca II IR lines 
are absorbed. However, estimating the non-LTE level popula- 
tions and the exact contribution of other lines to the observed 
profile needs detailed modelling of the SN spectrum. The high- 
est velocity of the three SNe as measured from the blue edges 
of their Ca II H&K profiles is shown by SN 1984A (38 000 
km s~') at -7 days (see Table 4). 

5.2. High-velocity Si II 

In Fig. 8, the Si II '6150 A' profile is shown for SNe 200 lei, 
1990N, and 1984A. Again, the data were normalised and con- 
verted to velocities w.r.t. the bluest component of the Si II pro- 
file of 6347 A (see Table 4). In Fig. 8e, the Si II profile of 
SN 200 lei at -9 days is compared with the Si II profiles of SNe 
1984A and 1990N at -7 and -14 days, respectively. At these 
epochs both SNe 200 lei and 1990N show a flat-bottomed Si II 
profile. To our knowledge such flat-bottomed Si II lines have 
not been seen in any other SNe la so far. 

The pre-maximum light evolution of the Si II profile is il- 
lustrated in Fig. 8f for SNe 200 lei and 1990N. For the latter 
epochs we used spectra presented by Wang et al. (2003) and 
Leibundgut et al. ( 1 99 1 ), respectively. In SN 200 1 el the flat bot- 
tom of the line disappeared by the next earliest epoch of obser- 
vation, 4 days before the maximum light. The same behaviour 
is also apparent in SN 1990N where the shape of the line has 
changed by the -7 day epoch. In Fig. 5a of Fransson (1984) 
such a flat-bottomed profile, and its evolution are presented as 
a result of calculations assuming that the velocity in the SN 
ejecta is proportional to its radius as expected soon after the SN 
shock outbreak. The flat-bottomed line shape is produced by 
pure scattering coming from a thin region compared to the size 
of the SN photosphere. As the scattering region widens when 
the SN 'photosphere' receeds the flat-bottomed line shape dis- 
appears, becoming an ordinary P-Cygni profile (cf., the change 
from the R2 = 1.2 case to the R2 = 1.5 case in Fig. 5a of 
Fransson 1984). 

The fact that the Ca II lines extend to considerably higher 
velocities than the Si II line (see Table 4) explains the less pro- 
nounced flat section of the Ca II lines. The lower maximum 
velocity of the Si II line compared to the Ca II lines is likely to 
be a result of the larger population of the atomic levels of the 
Ca II lines compared to the excited Si II line, and does not need 
to indicate different abundance distributions. We note that the 
flat-bottomed Si II profile in SN 1990N has also been explained 
by a blend of the Si II line with C II 6580 A line forming 
in a carbon-rich HV shell detached from the SN photosphere 
(Fisher et al. 1996), but see Mazzali (2001). Our explanation 
for the flat-bottomed Si II profile therefore gives less support 
to an enhanced HV carbon shell. 

6. Discussion and Summary 

We have looked for narrow hydrogen and helium lines in our 
early high-resolution optical spectra of the normal type la su- 
pernova SN 200 lei. No such lines were detected, and we de- 



Table 4. Maximum velocities for SNe 1984A, 1990N, and 
200 lei as indicated by the blue edges of the Ca II H&K and 
IR triplet and Si II '6150 A' profiles (see Fig. 8). 



SN 


Epoch 


Vca II.H&K 


Vca n,m 


V si n 




(days) 


(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


1984 A 


-T 


38 000 




26 000 


1990N 


-14* 


36 000 


-28 000 


26 000 


200 lei 


-9 


34 000 


28 000 


23 000 



" Date of observation 1984 Jan 10 (Wegner & McMahan (1987). 
The epoch adopting Jan 17 for the B-band maximum light (Barbon 
et al. 1989). 

b Date of observation 1990 June 26.2 UT (Leibundgut et al. 1991). 



rived upper limits for the luminosities of Ha, He I /15876 and 
He II /14686. These luminosities were then compared with 
the photoionisation models of Lundqvist et al. (2006). From 
the limits for Ha we conclude that the mass loss rate from 
the progenitor system of SN 2001el was less than < 9 x 
10" 6 M yr _1 for a spherically symmetric wind with a velocity 
of 10 km s _1 or less than < 5 x 10 5 M yr _1 for a velocity 
of 50 km s _1 . These estimates are only sensitive to gas out- 
side a few xlO 15 cm from the SN explosion site. The gas inside 
this radius would have been swept up by the SN ejecta prior to 
our observations. However, these are the best Ha based upper 
limits obtained for a SN la so far. Our results, together with 
the previous observational work on SN 1994D (Cumming et 
al. 1996) and SN 2000cx (Lundqvist et al. 2006) do not favour 
a symbiotic star in the upper mass loss rate regime (so called 
Mira type systems) from being the likely progenitor scenario 
for these SNe. This is in accordance with the models of Hachisu 
et al. (1999a, 1999b), which do not indicate a mass loss rate 
of the binary companion in the symbiotic scenario higher than 
~10~ 6 M yr _1 . Stronger mass loss would initiate a powerful 
(although dilute) wind from the white dwarf that would clear 
the surroundings of the white dwarf, and that may even strip off 
some of the envelope of the companion. The white dwarf wind 
with its possible stripping effect, in combination with the or- 
bital motion of the stars, is likely to create an asymmetric CSM. 
Effects of a more disk-like structure of the denser parts of the 
CSM, and how this connects to our spherically symmetric mod- 
els, are discussed by Cumming et al. (1996). In general, lower 
values of M/ V w ind should be possible to trace in an asymmetric 
scenario, although uncertainties due the inclination angle and 
the flatness of the dense part of the CSM would also be intro- 
duced. However, to push the upper limits down to the ~ 10~ 6 
M yr _1 level a much more nearby (D ~ 3 Mpc) SN la needs 
to be observed even earlier (~ 15 days before the maximum 
light) than our first epoch observation of SN 200 lei (Mattila 
et al. in prep.). This will only be feasible by making use of 
target-of-opportunity (ToO) time offered by observatories such 
as ESO and Subaru. Similar, or even lower limits on M/V w j n d 
can be obtained from radio and X-ray observations (e.g., Eck et 
al. 2002; Lundqvist et al. 2006), but such observations cannot 
disentangle the elemental composition of the gas in the same 
direct way as optical observations can do. 




Fig. 8. (a) Comparison between the Ca II IR triplet (solid line) and the Ca II H&K velocity profiles (dashed line) at -9 days, 
(b) Spectral evolution of the Ca II IR triplet from -9 days (solid line) to -4 days (dashed line), (c) The Ca II H&K profiles in 
SNe 2001el at -9 days (solid line), 1984A at -7 days (dotted line), and SN 1990N at -14 days (dashed line), (d) Comparison 
between the Ca II IR triplet profiles of SNe 2001el (solid line) and 1990N (dashed line). Note that due to the limited wavelength 
coverage in SN 1990N spectrum the normalisation of its Ca II IR profile is unreliable in the red part, (e) Comparisons between 
the Si II '6150 A' profiles in SNe 2001el at -9 days (solid line), 1984A at -7 days (dotted line), and SN 1990N at -14 days 
(dashed line), (f) Spectral evolution of the Si II '6150 A' profile in SNe 2001el (solid and dashed lines) and 1990N (dotted and 
dotted-dashed lines). The flux calibrated (f A ) spectral features have been normalised and converted to velocity space w.r.t. to the 
rest wavelengths assuming V rec = 1 180, 1010, and -261 km s -1 for SNe 200 lei, 1990N, and 1984A, respectively. For this we used 
the rest wavelengths of the bluest components of the Ca II triplet/doublet lines viz. 8498 A and 3934 A, and the Si II '6150 A' 
profile of 6347 A. 



SN 200 lei was revisited using VLT/FORS in the nebular 
phase, ~400 days after the maximum light. We modelled the 
late time spectrum to derive an upper limit of ~0.03 M Q for 
solar abundance material present at velocities lower than 1000 
km s _1 within the supernova explosion site. This is similar to 
our results for SNe 1998u and 2000cx (Lundqvist et al. 2006). 
Comparing this limit to the numerical simulations of Marietta 
et al. (2000) indicates that symbiotic systems with a subgiant, 
red giant or a main-sequence secondary star at a small binary 
separation are not likely progenitor scenarios for these SNe (see 
also below). These results demonstrate that the combination of 
very early and late time spectroscopy can be a powerful tool to 
probe SN la progenitors. 

Our -9 day spectrum shows that both the Ca II IR triplet 
and the H&K lines are present at very high velocities of up 
to -28 000 km s _1 and -34 000 km s _ \ respectively. High- 



velocity Ca II in SNe la is, however, not unique for SN 2001el. 
Here we have highlighted this for SNe 1984A and 1990N, 
which actually showed ~ 2000 - 4000 km s _1 higher velocities 
for the Ca II H&K as well as the Si II lines at a similar epoch. 
As a matter of fact, high-velocity Ca II has now been seen in 
all SNe la (Mazzali et al. 2005b) which have been studied early 
enough. The presence of the lines confirm the assumption in 
Lundqvist et al. (2006) that at least some SNe la launch their 
outermost ejecta at velocities in excess of c/ 10. This means that 
they sweep up any CSM at a much faster rate than was assumed 
in the earlier models of Cumming et al. (1996) (see Sect. 4.2). 

A fundamental question is of course the origin of the HV 
Ca II H&K lines. One possibility mentioned in Sect. 5 is 
that they form in material stripped off from a non-degenerate 
companion. From our FORS observation at 400 days post- 
maximum we were able to put an upper limit of ~0.03 M for 
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solar abundance material with a low velocity (~1000 km s _1 ) 
evaporated from the companion star. According to the mod- 
els of Marietta et al. (2000) only a small fraction of the gas 
lost from the companion, <1/1000, can have a high velocity 
of >15 000 km s _1 . The mass of the solar abundance HV gas 
(originating from the companion star) available to give rise to 
the HV Ca II lines would therefore be less than ~3 x 10~ 5 M 
which is a factor of ~3000 less than the mass of solar abun- 
dance gas estimated by Wang et al. (2003) to be required to 
produce the HV Ca II IR triplet feature seen in SN 200 lei. Our 
observations therefore do not support that the HV Ca II lines 
would arise in gas stripped from the companion star, which 
is consistent with the findings of Thomas et al. (2004) for 
SN 2000cx. 

We also note that there is no H/3 absorption apparent in our 
spectrum (see Fig.0 at the velocity of the HV Ca II features. 
A very tentative identification of such a HV H/3 line was made 
by Branch et al. (2004) in SN 2000cx spectra. The existence of 
a HV H/3 line would have indicated a non-degenerate compan- 
ion star origin for the HV matter, which is thus not supported 
by our observations. A very low hydrogen mass upper limit of 
~3 x 10~ 4 M was presented by Delia Valle et al. (1996) for 
SN 1990M based on a non-detection of atmospheric Ha ab- 
sorption in their spectrum near the maximum light. However, 
this was obtained by using a simplified model. Ha is usually in 
net emission in supernovae, and thus the absorption is weaker 
than would be predicted by an analysis using resonant scat- 
tering. The emission falls into the Si II absorption trough and 
so could be hidden in SNe la (Thomas et al., 2004; Branch et 
al. 2004). This effect is somewhat less for H/3 and Hy. Lentz 
et al. (2002) performed detailed spectral synthesis calculations 
of mixed solar material into a W7 model and found that even 
relatively large amounts of hydrogen (~ 10~ 3 M Q ) would be 
hard to detect at the earliest times. Furthermore, early time ob- 
servations can only detect hydrogen present in the outer parts 
of the supernova ejecta, or in a detached wind (see Sect. 4.2). 
However, as the simulations of Marietta et al. (2000) predict 
only a very small high-velocity tail for the hydrogen rich mate- 
rial, early observations seem futile to reveal any stripped mate- 
rial from the companion. The models of Marietta et al. (2000), 
in combination with our results in Sect. 4.3 and in Lundqvist 
et al. (2006), show that stripped, or rather evaporated material 
from a companion, should instead be looked for at late times. 
Our limits on the late emission are actually so tight that most 
of the scenarios in Marietta et al. (2000) appear to be ruled out. 
This, in combination with the properties of the presumed bi- 
nary companion of the Tycho progenitor (Ruiz-Lapuente et al. 
2004), should be used to find a viable solution to the single- 
degenerate progenitor scenario. 

Another possibility for the HV Ca II lines mentioned in 
Sect. 5 is that they arise in gas compressed due to circum- 
stellar interaction (Gerardy et al. 2004). The scenario with a 
dense continuous wind outside a fewxlO 15 cm is not supported 
by our observations and models of narrow circumstellar emis- 
sion lines. Furthermore, the lack of interaction signatures in the 
bolometric light curves (e.g. Gerardy el al. 2004) and the lack 
of radio emission in normal SNe la (e.g., Eck et al. 2002), do 
not support such a scenario either. The model with a circum- 



stellar disk could be more likely, although we note that such 
a scenario should produce a range of velocities for the Ca II 
absorption depending on the viewing angle, unless the disk is 
thick. More self-consistent modelling of the circumstellar in- 
teraction is needed to see whether a circumstellar disk model 
can explain the HV Ca II features seen in SNe la. 

A more natural explanation for the HV Ca II lines and 
the HV Si II lines in SN 200 lei and other SNe la would ap- 
pear to be absorption in the outermost parts of the SN ejecta 
which extend to high enough velocities. We discussed this qual- 
itatively in Sect. 5.2. Maximum velocities of up to ~26 000 
km s _1 and ~38 000 km s _1 were observed for the blue edges 
of the Si II and Ca II H&K profiles, respectively, in the -7 
day spectrum of SN 1984A. Lentz et al. (2001) modelled the 
early spectra of SN 1984A and found that the large blueshifts 
could be produced by a higher density in the outer layers of 
the ejecta. Furthermore, Hatano et al. (2000) showed that even 
at 10 days past maximum light, i.e., much later than the times 
studied here, there exists a class of 'fast' supernovae typified 
by SN 1984A. More recently Benetti et al. (2005) divided 'nor- 
mal' SNe la into groups showing high (HVG) and low (LVG) 
temporal velocity gradients for the Si II profile after the max- 
imum light. Both SNe 200 lei and 1990N belong to the LVG 
group whereas SN 1984A belongs to the HVG SNe which 
also tend to show the highest photospheric velocities both be- 
fore and after the maximum light. It is important to understand 
whether the density structure variations of the HV SNe la are 
due to initial conditions, variations in the explosion mechanism 
or something else. How these HV SNe la are correlated with 
brightness at maximum light is therefore an important question 
for the use of SNe la as cosmological probes. 
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